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ABSTRACT 

Winds from short-period Earth and Neptune mass exoplanets, driven by high en¬ 
ergy radiation from a young star, may evaporate a significant fraction of a planet’s 
mass. If the momentum flux from the evaporative wind is not aligned with the 
planet/star axis, then it can exert a torque on the planet’s orbit. Using steady-state 
one-dimensional evaporative wind models we estimate this torque using a lag angle 
that depends on the product of the speed of the planet’s upper atmosphere and a flow 
timescale for the wind to reach its sonic radius. We also estimate the momentum flux 
from time-dependent one-dimensional hydrodynamical simulations. We find that only 
in a very narrow regime in planet radius, mass and stellar radiation flux is a wind 
capable of exerting a significant torque on the planet’s orbit. Similar to the Yarkovsky 
effect, the wind causes the planet to drift outward if atmospheric circulation is pro¬ 
grade (super-rotating) and in the opposite direction if the circulation is retrograde. 
A close-in super Earth mass planet that loses a large fraction of its mass in a wind 
could drift a few percent of its semi-major axis. While this change is small, it places 
constraints on the evolution of resonant pairs such as Kepler 36 b and c. 

Key -words: planets and satellites: atmospheres planets and satellites: physical 
evolution ultraviolet: planetary systems ultraviolet: stars X-rays: stars. 


1 INTRODUCTION 


Planets discovered by the Kepler Observatory 
llBatalha et al.|[2013l) exhibit a large range of d ensities rang¬ 
ing fr om less than 1 g cm~^ (e.g. Kepler 36c: ICarter et al.l 
I 2 OI 2 I and Kepler 79d, Ijontof-Hutter et al.l 201^ to gr eater 
than 8 g cm“^ (e.g. Kepler 10b, Batalha et ahlEoill) . For 
Neptune and Earth mass planets, the inferred gaseous 
envelope masses range from zero to tens of percent of the 
total planetary mass, with this fraction bein g dependent 
on w hether a rock/iron or ice core is adopted (iHowe et al.l 
[20 i 3). Either close-in planets are formed with a wide range 
of core to envelope mass ratios, or their compositions evolve 
after they have accreted gas (and after evaporation of the 
circumstellar disk). Evaporative winds driven by stellar 
UV and X-ray radiation could account for the loss of a 
significant frac tion of a Neptune-sized planet’s gaseou s 
envelope mass (iLopez fc FortneJl2013l : lOwen fc Wull2013l) . 
Impacts or collisions between bodies that are more frequent 


at small semi-m ajor axis can also affect the planetary 
composition (e.g., fouillen et ahllioU) . Among planet pairs 
in multiple-planet systems, a larger fraction (~ 60%) of 
systems have the smaller radius planet residing closer to 
the star dCiardi et al.ll2013l). Along with the clos er in planet 
generally being denser i Wu fc LithwickI l2013l) . As more 
mass can be evaporated clos er to a star, evapor ative winds 
could account for this trend llOwen fc Wull2013) . 


In addition, evidence for mass loss through atmo¬ 
spheric evaporation h ave been observed for the h ot Jupiters 
around HD 209458 (fW dal-Madiar et al~1 l2003l) and HD 
189733 (iLecavelier des Etangs et al. 2010l : Bourrier et al.l 
I2OI3I) , with loss rates estimated to be on the order of 
10i° gs-k 


An evaporative wind from an exoplanet is often as¬ 
sumed to be isotropic. However, an evaporative wind car¬ 
ries momentum and could exert a torque on the planet 
affecting both its orbit and spin. Here we first summa- 
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2 Teyssandier et al. 


rize related physical phenomena. IVeras et al.l (|201.ll l showed 
that an anisotropic stellar wind can cause orbital inclina¬ 
tion changes. The Yarkovsky effect occurs when uneven re¬ 
emission of absorbed stella r radiation causes drifts in as¬ 
teroid semi-major axes Isee lBottke et al]l200fil l. It has also 
been shown that sublimation from a c omet can exert a reac¬ 
tive torque on the cometary nucle us (iG utierrez et al.l I 2 OO 2 I : 
iNeishtadt et al.ll2003l 'l. iBoue et ^ ll2012|j proposed that an 
anisotropy in an evaporative wind from a close-in Jupiter- 
or Neptune-m ass planet c an cause the plane t’s orbit to mi¬ 
grate (see also llori a l2012l ~). iBoue et ^ (120121 ') suggested that 
migration distances inferred from distributions of planetary 
mass and orbits could be used to constrain the launch angle 
of the evaporative wind. In this paper we attempt to esti¬ 
mate the launch angle of the wind and the magnitude of 
the resulting reactive torque on the planet based on hydro- 
dynamical models. As evaporation is primarily significant 
for low core mass planets experiencing high UV and X-ray 
fluxes llLopez fc Fortne^1l2013l : lOwen fc Wull20l3 i from the 
star, we focus here on winds driven from warm Neptunes 
and super-Earth like planets within 0.2 AU of a star. 

As the cooling timescale in an ionosphere is expected 
to be shorter than a day, an evaporative wind is mainly 
expected to be driven from the day side of a planet. If the 
planet is not tidally locked or if the atmosphere is circulating 
then the upper atmosphere would experience day/night vari¬ 
ations in stellar illumination. A delay between dawn and the 
time it takes for the wind to develop could cause an asym¬ 
metry in the wind (see Figure [T]). Here we aim to estimate 
the torque on a planet due to this possible asymmetry. 

Planets within 0.2 AU are expected to be tidally locked 
and at low obliquity (e.g., iLubow et al.l[l9^ : ICunha et al.l 
|2014) . Spin-down sh ould have taken pla ce in less than ~ 10^ 
years (see Table 1 bv ICunha et al.|[20l4 ). and well before the 
time when most of the evaporation takes place, at an age 
~ 10^ — 10® years (lOwen fc Jacksonll2012ll . However, some 
planets could now be in a spin synchronous state (with rota¬ 
tion rate related to the angular rota tion rate at pericenter , 
iHutI llOSd 'l or spin-orbit resonance (iRodriguez et al.l I 2 OI 2 I ') 
like Mercury. Even though the tidal spin-down timescale 
is short, the timescale for tidal circularization of the or¬ 
bit i s relatively long (e.g., iFord &: Rasid l2008l : I Lee et akl 
l2013l 'l and planets could have been at moderate (or their 
current) eccentricities when atmospheric evaporation took 
place. About a quarter of the planets in multiple-planet Ke¬ 
pler systems with transit-timing variations have high ec cen- 
tricities, in the range e=0.1-0.4 (IWu fc Lithwickll2013lj . In 
exoplanetary systems containing at least two planets, tran¬ 
sit timing variations are consistent with an eccentricity dis- 
ersion of O.OlSln nnj whe n fit with a Rayleigh distribution 


Hadden fc Lithwick 


2OI4II . 


Even if a planet is tidally locked, the upper atmosphere 
can be rotating wi th advective or zonal flow speeds of order 
Vzonai ~ 1 km/s (Ishowman fc Polvanil I 2 OIII : iMenoull2012l : 
iRauscher fc Kemotonl 20141 ). If the advective flow sp eed is 
of order the sound speed ('e.g.. l^wan fc Agolll20lTbl 'l then 
we expect a wind speed oc T^q , where T^q is the equilibrium 
temperature set by radiation balance from the stellar light. 
As this is weakly dependent on the planet’s semi-major axis 
(Teq oc with ap the semi-major axis of the planet), 

for simplicity we can neglect the possible dependence of ad¬ 
vective flow speed on stellar flux. We also neglect the de¬ 


pendence on planet properties such as atmospheric opacity 
and surface gravity, as atmospheric models all predict simi¬ 
lar circulation speeds. Evaporative winds are driven from a 
tenuous region of a planetary atmosphere where ionization 
and heating occur, and this layer is higher than typically 
modelled to predict atmospheric circulation. Nevertheless, 
~1 km/s size-scale zonal fl ows are predicted in exoplanet 
ionospheres (see Figure 3 bv iKoskinen et ^l20ld j. 

For completeness, we note that, in addition to the ef¬ 
fects outlined above, magnetic fields on the planetary surface 
will lead to anisotropy. More specifically, with the expected 
mass loss rate from planets and surface field strengths of 
^ 1 Gauss, the outflo ws could be magnetica lly controlled 
(lOwen fc Adamj|2014l . see also lAdama ll201ll ~)~l. This issue 
must be addressed in future work. 

Both zonal circulation or a spin-orbit resonance would 
cause a planetary ionosphere to experience day/night varia¬ 
tions. In section^we estimate a torque due to an evaporative 
wind by considering the relationship between atmospheric 
motion and the timescale to launch a wind at dawn, as this 
is expected to dominate the rotation of gas in the upper 
atmosphere, unless the planet is in an unusual spin-orbit 
resonance. In section [3] we compute a lag angle for the wind 
momentum loss rate using one-dimensio nal hydrodynami- 
cal st eady-state evaporative wind models (I Owen fc JacksonI 
l2012|j . In section [4] we compute this angle by integrating 
one dimensional but time-dependent hydrodynamic mod¬ 
els. Lastly in [S] we investigate scenarios for the Kepler 36 
planetary system taking into account an evaporative wind 
induced torque. Variations in the planet’s orbital elements 
caused by a torque from an anisotropic evaporative wind, 
integrated over the orbit, are computed in the appendix. 


2 TORQUE ESTIMATE AND LAG ANGLE 

In this section we roughly estimate a lag angle and the 
torque exerted by the wind onto the planet. We first de¬ 
rive a launch timescale for the wind, and compare it to the 
circulation timescale in the upper atmosphere. 


2.1 Relevant timescales 


Stellar X-ray and EUV radiation heats ga s up to temper- 
atures T ^ K (ISpitzeil[l978l : [Shulll992l : [Ercolano et al.l 
I 2 OO 8 II corresponding to a sound speed of ceuv ~ 10 km/s. 
Above this temperature radiative cooling is so efficient that 
higher temperatures are rarely reached. As ceuv is similar 
to a planet’s escape velocity, a wind can be launched from 
the heated gas. An isothermal or Parker-type wind model 
(|ParkerjIl958l . 1 19651') gives an approximate hydrodynamical 
description for evaporative winds from close-in planets, i.e., 
planets for which evaporati ve winds can be impor tant (out 
to about 0.2 AU, see, e.g., lOwen fc JacksonI I 2 OI 2 II . Such a 
wind cannot be instantaneously launched at dawn. We can 
estimate a launch timescale for the wind from the time it 
takes a particle in the wind to travel from the planet’s sur¬ 
face, at a planetary radius Rp, to the transonic radius, Ra 
(where the flow becomes supersonic), and is no longer in 
causal contact with the planet. 


tflo 


dr 


Ra 


Ceuv 


( 1 ) 
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Figure 1. If the wind launch region is super-rotating with respect to the planetary surface, then the wind is launched after the upper 
atmosphere sees dawn. A comparison between a timescale to launch the wind and with the time required for the upper atmosphere to 
rotate gives an angular lag in the peak of the day-driven wind. Because of the lag, the momentum flux of the wind is offset compared 
to the star/planet axis and this gives a torque on the planet. If the upper atmosphere is super-rotating in the prograde direction with 
respect to the planetary surface, the torque increases the angular momentum of the planet and the planet’s orbit moves outward. 


where u is the wind velocity as a function of radius. The one¬ 
dimensional hydrodyna mic models that i nclud e ionization 
and heating by X-rays, I Owen fc JacksonI (l2012ll have Rs ~ 
2 — 3Rp (see their Figure 4) and ceuv ~ 5 km s“^, where 
Rp is the planet’s radius. Using these larger transonic radii 
we find: 


tflo 


Rs 


Ceuv 


10® s 


fRp\ 

(Rs/RA 

( Ceuv \ 

V^®y 

V 2 ) 

V 5 km s“® / 


( 2 ) 

The motion of the upper atmosphere wind, described by 
a zonal wind speed Vzonai, and the timescale to launch the 
wind together imply that the wind is not launched at dawn 
but slightly afterwards. It would take a similar time for the 
evaporative wind to die down after dusk. The delay implies 
that the average wind momentum is not pointing exactly 
toward the star along the sub-stellar point but is shifted by 
an angle Oiag from this direction (see Figure[T]). We estimate 
this angle to be 


9 la 


tflo 


t zonal 


where 


tz 


Rp 


= 637s 


f Vzonai 


Vlkm s-W 

U®y 


(3) 


(4) 


'^zonal 

is the inverse of the atmosphere’s angular rotation rate. A 


comparison between equation @ and equation © suggests 
that the angle 9iag could be of order 1 and the wind could 
exert a significant torque on the planet. 


If tfio-w S> Rp/vzonai then we could consider the wind 
launched evenly from both day and night sides. If tfio-w <C 
Rp/Vzonai then the momentum flux from the wind (averaged 
over all regions on the planet surface) is in the radial direc¬ 
tion toward the star, and there is no torque exerted by the 
wind onto the planet. We have used a flow timescale to esti¬ 
mate the time it takes to launch the wind that is based on a 
steady state flow model. An alternative thermal timescale, 
tkh, that can be used to estimate a wind launch timescale is 
the time it takes to replenish the kinetic energy in the flow 
within the sonic radius from the absorbed EUV and X-ray 
energy flux. Depending upon the ratio of the transonic to 
planet radius, and as we discuss below, a thermal timescale 
may more accurately describe the time it takes to launch 
the wind. We detail this possibility in section [3] We have 
neglected the wind speed’s tangential velocity component in 
our estimate for 6iag as Vzonai is significantly smaller than 
the escape velocity from the planet surface. 
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2.2 Orbital and spin torques exerted by the wind 

The orbital torque (from the wind) can be roughly estimated 
from the lag angle 9iag as 

Till ~ ^GlagVuidj)/‘I ( 5 ) 

where M is the wind mass outflow rate and Vm ^ Ceuv is 
a mean wind velocity near the transonic region. The rough 
factor of 2 comes from integrating a radial wind over a hemi¬ 
sphere with an angular offset. 

The sign of the lag is set by the direction of the planet’s 
spin or atmospheric flow. If the planet’s spin or the at¬ 
mospheric circulation are prograde (super-rotating or east¬ 
ward) then the torque is positive and we expect the planet 
to move outwards. The sign is in dir ect analogy to the 
Yarkovsky effect fsee lBottke et aUl2006l 'l. with outward mi¬ 
gration in semi-major axis for prograde asteroid body rota¬ 
tion and vice versa for retrograde body rotation. Close-in ex¬ 
oplanet atmospheric circulation models usually predict pro¬ 
grade (super-rotati ng or eastward) equatorial upper atmo¬ 
spheric circulati on (ISho wman fc Polvanill201ll : lMenoij[2fll2l : 
iRauscher fc KemptonI 2014l b with the exception of the ret- 
ro grade (westward) slowly rota ting model for HD 209458b 
bv iRauscher fc Kemptor] ll2014ll . Although, we caution that 
the atmospheric rotation of low-mass close in exoplanets has 
not been studied and note in passing that the circulation di¬ 
rection in Neptune and Uranus is retrograde. 

Here we also note that the evaporative wind could affect 
the spin evolution of the planet. The torque on the planet 
itself we estimate as 

Ts ^ ^{Vzonal Rp^p) Ftp ( 6 ) 

where Vzonai is the atmosphere circulation velocity and we 
assume that the wind is launched at a radius that is not 
significantly different than Rp. Here Hp is the spin of the 
planet. Assuming that Vzonai > RpFlp and after time At, 
one gets a change in planet spin of 

MAt Vzonai 1 




Mp Rp 


(7) 


where the planet’s moment of inertia is Ip = OpMpRp. For 
short-period, tidally locked planets, this effect could act to 
desynchronise the planet’s spin, provided synchronisation 
occurred before the evaporation became important. 


2.3 Orbital evolution 

A strong evaporative wind is driven when the host star has 
a strong EUV a nd X-ray flux or durin g an age less than 
At ~ 10® years JOwen fc Jacksonll2012l '). During this time, 
a planet of mass Mp can experience a loss of mass AMp = 
MAt and a change in semi-major axis Adp. The change in 
the planet’s orbital angular momentum is 

AL « Till At ~ MpVcAdp/2, (8) 


where Mp is the planet mass and Vc the planet’s orbital 
velocity (assuming a nearly circular orbit), which scales as 


Vc = 94km s ^ 


/M* \ 2 / ap 
\Mq) VO.IAU/ 


Using equation m for the torque, we estimate 


^^Qjp 

dp 


MAt Vu 

Mo Vc 


■6la 


A]\Ip Viii 

Mo Vc 


lag 


(9) 


( 10 ) 


The factor Vmjvc limits the size-scale of the possible migra¬ 
tion distance in semi-major axis. 

Many o f the Kepler planet h ost stars have masses lower 
than IMq (iBatalha et al.l l2013l l and for these lower mass 
stars the ratio Vmlvc would be higher. Previous studies 
jOwen fc WuII20i 3 : iLopez fc Fortne^i 2013l 'l have found that 
close-in low core mass planets can loose a significant fraction 
of their mass. Taking AMp/Mp ~ 1 and 9iag ~ 1, Equation 
(fTOt gives 


Adp ^^AMp^ ( Vm \ (M,\'^ ( dp 

dp ^ Mp '“nio kms-iHA70y VO.IAU/' ’ 

( 11 ) 

A planet that loses a significant fraction of its mass could 
move a significant fraction (say 10%) of its semi-major axis. 
10 % exceeds the distance to resonance exhibited by resonant 
pairs l^brvckv et al.l ll2014h and is dynamically significant, 

J iarticularly for clo se resonant pairs such as Kepler 36b,c 
Carter et aLllioi^ i. 

In appendix]^ we derive more general equations for the 
secular evolution of semi-major axis and eccentricity for a 
planet whose orbit is perturbed by an evaporative wind. 


3 LAG ANGLE ESTIMATED FROM 

STEADY-STATE MODELS 

In the previous section we estimated the torque on a planet 
due to evaporative wind loss using a crudely approximated 
launch timescale for the wind. In this section we use the 
steady-state hydrodyn amic models, including r adiative cool¬ 
ing and ionization, bv lOwen fc JacksonI (l2012l l to better es¬ 
timate the lag angle 9iag that sets the torque from the wind. 

As discussed above we take the flow timescale to be the 
time to establish an evaporative flow. It is the time the flow 
takes to reach the sonic point and is thus found by: 


tfi 


OW 



dr 

u 


( 12 ) 


wher e u is the velocity field of the flow. The flow solutions 
from lOwen fc JacksonI (l2012l i allow us to directly perform 
this calculations without making further assumptions. A 
comparison between this flow timescale and the zonal wind 
timescale using a zonal wind velocity of 1 km s“^ (Equa¬ 
tion [4]) is shown in the left-hand panel of Figure (2] In line 
with our discussion in Section 2 the flow timescale is of order 
the zonal flow timescale for Neptune like planets. 

One further consideration is that the flow receives 
enough energy to launch the wind, we call this timescale the 
“Kelvin-Helmhotlz timescale” of the flow (following Owen 
& Adams, 2014). We can also estimate the time it takes to 
heat the wind from a ratio of kinetic energy (integrated out 
to the sonic radius, ~ psCs47rR®/3) and the energy absorbed 
per unit time in UV and X-ray photons, ~ ixR^Fuv where 
Ps is the density in the wind at Rs and Fuv is the flux in 
X-ray and EUV photons from the star. In analogy to the 
Kelvin-Helmholtz timescale we write 


tkh 


A-KT^Arpy? 
TvRpFuv ’ 


(13) 


where the numerator represents the integral of the energy 
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density in the flow over the volume within the sonic radius. 
We show the calculated “Kelvin-Helmhotz timescales” com¬ 
pared to the zonal flow timescale in the right panel of Fig¬ 
ure [2l We And that the flow timescale dominates the launch 
timescale when the potential is deep (as the suppression of 
the launch velocity is exponentially sensitive to planet mass) 
and the flow is not close to “energy-limited” (see discussion 
in Owen & Jackson 2012). However, when the potential is 
shallow, the sonic point is very close to the planet’s sur¬ 
face and the flow timescale is essentially zero. In this case 
the timescale to launch the wind is controlled by the energy 
considerations. 

In reality both timescales matter. Therefore we define 
a “launch timescale” as the root-mean square of the tow 
timescales, and we compare it to the zonal flow timescale. 
This comparison is shown in Figure [S] It shows that there 
is a reasonably wide band of planet parameters where the 
launch timescale compared to the zonal timescale is order 
unity. We compare it to the existing sample of exoplanets. 
We use the database from the Open Exoplanet CatalogutQ 
n einll2012i l. for all planets with orbital period lower than 
100 days. A small fraction of the super-Earth to Netpune- 
like planets are in the range where the timescale ratio (and 
hence our estimated lag angle) is of order unity. 


4 MOMENTUM LOSS RATE ESTIMATED 
FROM ONE-DIMENSIONAL TIME 
DEPENDENT HYDRODYNAMICAL 
MODELS 

To improve on the estimates of section [3] we carry out some 
one dimensional but time-dependent hydrodynamic models. 
From these we directly estimate the lag angle and the vector 
direction of the momentum carried by the evaporative wind. 

4.1 Method 

We conduct a set of simulations using ZEUS 
dStone &: NormM ] Il992f) . The code we us includes 
modifications for UV ionization described in detail by Owen 
& Adams (2014) as well as time-dependant heatin g and 
cooling described in detail bv lOwen fc Alvared (l2015l i. The 
numerical method takes into account ionization by a range 
of EUV photons, and as such includ e s pho ton-hardening 
using the method of iMellema et al.l ll2006l l. assuming a 
pure Hydrogen atmosphere, where ray-tracing is performed 
in a causal manner starting at the edge of the grid and 
working towards the planet. Heating is dominated by 
photo-ionization and cooling from Lyman-tr and recom¬ 
binations. 3D simulations of this problem are beyond the 
scope of this work and therefore we chose to setup a simple 
illustrative problem, that combined with the insight from 
our discussion above, allows us to make inferences about 
the role of evaporation in driving a planetary torque. 

Our setup is as follows: we consider a small patch 
of the planet rotating with azimuthal velocity v,p such as 
"Uzonai = cos 8 where 6 is the latitude of the patch. For 
each patch we solve the U-ID (including the evolution of 

^ http://www.openexoplanetcatalogue.com/ 


v^) in spherical co-ordinates. We note that this implicitly 
assumes that the flow has spherical divergence along the 
streamline and we neglect any pressure forces perpendicu¬ 
lar to the radial direction. This is obviously a simplification; 
however for a thermally driven outflow the pressure forces in 
the radial direction must dominate. During the calculation 
we then modulate the input UV flux as a function of time 
to mimic the patch rotating from the day to night side and 
wait to a quasi-repetitions flow is established, after several 
rotation periods of the patch. Once we have simulated dif¬ 
ferent patches at different latitudes we can then “stitch” the 
flow patterns together to gain an estimate of what the 3D 
flow structure is like. In particular, we can then explicitly 
estimate what the net torque on the planet is. As we have 
stitched together several patches one loses the constrain the 
momentum flux will be conserved in a global 3D sense; and 
as such we must be cautious in interpreting our results. 

The grid is radially divided in 128 cells, and the time- 
step is such that it takes about 100 iterations for a patch 
to go around the planet. In total there are 61 patches going 
from south pole to north pole. The outer boundary of the 
radial grid is such that all the flow becomes supersonic inside 
the grid. 

4.2 Results 

First, we study a Neptune-like planet, with a radius of 2.4 x 
10® cm and a mass of 1.03 x 10®® g. It is irradiated with a 
flux of 2 X 10® erg s“^ cm~®. This corresponds to a planet 
at 0.013 AU from a star with a X-Ray and UV luminosity 
of 10®® erg s“®, i.e. a young solar-type star. In what follows, 
we also note p, P and u the density, pressure, and velocity 
fields of the flow. 

In figure |4] we show the sonic point for different lati¬ 
tudes and longitudes, on the day-side only. As expected, the 
sonic point is close to the planetary surface around the sub- 
stellar point. However, the minimum of the curve appears 
to be offset from the substellar point by a few degrees. We 
estimate that, at the equator, the minimum sonic radius is 
offset from the sub-stellar point by 16° westward. 

From this set of simulation, we can also extract the 
total mass loss, from our combined stitched together 3D flow 
solution, integrated over the longitudes 0 and latitudes 8, at 
an altitude r (e.g., the sonic radius). It reads: 

M= f [ sm8Add4>p{r^8,(j))u(r,8,(f>)r^. (14) 

Je=--n/2 J,j>=0 

In this case we And it to be 1.3 x 10®® g s ~®, similar to the 
ID sp herically symmetric mass-loss rate of lOwen fc JacksonI 
( 1201 ^ . 

Another relevant quantity is the the momentum loss 
rate of the wind, P, whose components i = x, y, z are 
given by: 

Pi = — [ puiUjdA^ — [ PdAi — f dVpditpg, (15) 
J dV JdV Jv 

where ijig is the gravitational held of the planet. P is a vector, 
which we give in the following reference frame: the origin is 
at the substellar point, the x direction lies along the line 
that goes from the substellar point the the star, the y axis 
is perpendicular to the x axis along the equator in the east 
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Figure 2. Left: as given by eq. ((Tg} for a n X-ray + EUV stellar luminosity of 10^® erg/s and a semi-major axis of 0.1 AU, 

using the results of numerical simulations presented in lOwen &: Adamd ||2014| '). Right: same for tfc/i/tjonaZ given by equation II13II . In 
both plots, the star represents Neptune, while the circle represents a pre-evaporation Kepler 36b. 
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Figure 3. Root mean square of tfiom/tzonal as given by eq. I I12II 
for an X-ray -f EUV stellar luminosity of 10^® erg/s and a semi¬ 
major axis o f 0.1 AU, using the res ults of numerical simulations 
presented in lOwen fc AdamsI ll2014l ). and tkh/tzonal as given by 
equation [13] As a guide, we also plot a sample of exoplanets with 
orbital periods less than 100 days. In addition, the star represents 
Neptune, while the circle represents a pre-evaporation Kepler 36b. 


direction, and the 2 is perpendicular to the equator in the 
north direction. The problem is symmetric with regard to 
the planet’s equator, so that the 2 contribution cancels out 
once integrated over the whole planetary surface. 

We estimate the lag angle of the wind 9iag through the 
relation 

tan = Py IPx- (16) 

The integration is made in spherical coordinates, over our 
3D stitched together flow solution, at a constant distance 
from the planet surface. 

In the present simulations we find that the value of 
diag depends on the distance at which the integration is 
conducted, as shown on figure [S] This is a consequence of 
the fact that our stitched flow solution does not conserve 
momentum explicitly in the angular directions as discussed 
above. In any case, we find that the mass loss is strongly 
anisotropic, with a lag angle of several tens of degrees. 

Our one dimensional time-dependent model has there¬ 
fore allowed us to show that a short-period, strongly irra¬ 
diated, Neptune-like planet, looses mass at rate similar to 
other studies. In addition the momentum flux indicates that 
mass loss is anisotropic, and occurs at an angle that is offset 
from the substellar point by a few tenths of radians. Ac¬ 
cording to our estimate in equation (nnj, this would shift 
the semi-major axis of the planet by about 0.3% after 100 
Myr. In appendix El we also give a more general formula 
for the rate of change in semi-major axis due to the force 
exerted by the wind, which we can extract directly from the 
momentum flux. Using eq. (TAill we find that over 100 Myr, 
the hot Neptune that we considered above would have mi¬ 
grated outward by about 1% from the location at which we 
estimated the momentum flux (0.013 AU). 


© 0000 RAS, MNRAS 000, 000-000 


zonal 


























Torque on a planet from an evaporative wind 7 



0 71/4 nil 371/4 71 

longitude (rad) 


w 

3 

Ti 

CO 

DC 

-I—• 

CD 

C 

_C0 

Q. 


c 

o 

D_ 

o 

'c 

o 

CO 


Figure 4. 2D map of sonic radius divided by planetary radius for a Netpune-like planet, as a function of longitudes (day-side only, 0 is 
dawn) and latitudes (from south pole to north pole, 0 is at equator). The minimum sonic radius is offset from the substellar point by an 
angle of approximatively 16° westward. 



Figure 5. Lag angle 0 as a function of the radial distance r from 
the planetary surface, for a Neptune-like planet. 


5 APPLICATION TO KEPLER 36 

Kepler 36b has a density pi, = 7.46 g cm“^ whereas Kepler 
36c has Pc = 0.89g cm“®. Kepler 36b is the lighter denser 
planet, suggesting that it could have lost much of its mass in 
a wind. The Kepler 36 exoplanet system is thought to have 
obtained such a dichotomy in densities du e to evaporation 
llOwen fc Wull201^ : iLopez fc Fortnevll2013l ) . 

5.1 Torque due to hydrodynamical escape 

If the torque on the planet was positive during evaporative 
wind loss then the planet would have migrated outward. 


pushing it closer to Kepler 36c. This would have pushed the 
system deeper into the 7;6 resonance rather than pulling it 
slightly out of resonance. 

Currently, the system is just slightly outside of reso¬ 
nance with the separation be tween the two bo dies larger 
than that of exact resonance llCarter et al.l[201^ 1. The two 
planets in this system could not have migrated much fur- 
ther apart as tha t would have made the system unstable 
dDeck et al.ll2012l ). 

iLopez fc FortnevI (l2013l ) have found that the peculiar 
architecture of the Kepler 36 system can be reproduced by 
evaporation of a pair of planets with cores of 4.45 (in¬ 
ner planet) 7.34 and (outer planet). Both cores are sur¬ 
rounded by a gaseous atmosphere, which initially contains 
about 23% of the total mass of the planet in the form of 
hydrogen and helium. Because of the close proximity of the 
pair of planets and the star, the atmospheres of the form¬ 
ers undergo hydrodynamical escape, mainly during the first 
10® yr, when the stellar flux is stronger. Since the core of 
the inner planet is less massive, it is unable to keep its at¬ 
mosphere and only the dense rocky core remains, which we 
observe today. The outer planet manages to retain an ex¬ 
tended gaseous envelope, resulting in an observed low mean 
density. The two planets of the Kepler 36 system are also 
very close to a 7:6 mean motion resonance. In this paper 
we study the possibility for the inner planet to have been 
slightly shifted out of the resonance under the action of the 
torque exerted by its own evaporating atmosphere. 

We conduct hydrodynamical simulations similar to the 
ones presented above, for a planet consisting of a rocky 
core of 4.45 and an atmosphere whose mass is 22% 
of the core mass, and extends to 10 R^, in agreement to 
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the model bv lLopez fc FortnevI (|2013h . It receives a flux of 
2.2 X 10"^ erg s“^ cm“^ at 0.12 AU. We the perform a series 
of 1+lD simulations to map out the flow structure using the 
method described in Section [H with the same resolution. 

The sonic radius as a function of longitudes, at the equa¬ 
tor, is presented in £gure[6l and includes both day and night 
sides. The minimum sonic point is largely offset from the 
substellar point, in the East direction. On the night side, 
the atmosphere does not entirely cool down and a wind still 
exists because of the temperature at the base of the atmo¬ 
sphere, which is set to 1000 K in our simulations. 

On figure [3 we observe that the lag angle is minimum 
(and negative) at intermediate radii (r ~ 7 rpianet)- In our 
simulations we remark that on the days-side, the wind is 
rapidly launched and decreases in density at large radii and 
never reaches a quasi steady state. On the night-side, there 
still exists a transonic point at which the wind is launched 
(see figure [6ll. This wind is in quasi steady-state and isother¬ 
mal, and has a larger density at large radii than the day-side 
wind. This can possibly explain the lag angle profile of figure 
[3 This is an example of a case where multi-dimensional sim¬ 
ulation are required in order to fully understand the prob¬ 
lem. 

We find the mass loss to be 3.7 x 10^^ g s“^. On figure[3 
we show that the lag angle strongly depends on the position 
at which it is evaluated, but does not exceed 8°. This lag 
angle is smaller than the one we found for a Neptune-like 
planet, which conhrms the simple estimate we made in sec¬ 
tion |3 (see left panel of £gure[3l- Figure |T] clearly indicates 
that for the Neptune-like planet, the sonic point was located 
in a small range from 1.35 to 1.5 planet radius. In the case 
of Kepler 36, the sonic point spreads a much wider range of 
values. If we ignore the contribution from the polar regions, 
most of the region with r/vp > 8 is transonic. In this region, 
the lag angle increases with r, just as in the Neptune-like 
case. From equation (nni) we estimate that the resulting drift 
in semi-major axis would be, at best, of 0.4% after 100 Myr. 
This is fairly consistent with equation (IX4l) . which predicts 
an inward drift of 0.2% after 100 Myr when the torque is 
maximum. 

We have also conducted simulations for a pre- 
evaporation Kepler 36c planet, using again the model by 
iLopez &: Fortne^ (l2013l j. We find that the lag angle increases 
from 15° to 45° degrees, as we increase the radial location at 
which the momentum flux is evaluated. The mass loss rate 
is 1.3 X 10^^ g s“^. Using equation (nsj we estimate that the 
resulting drift in semi-major axis would be of about 0.2% 
after 100 Myr. 


5.2 Effect on the resonant pair 

The period ratio of planets b and c in the Kepler 36 system 
is Pc/Pb = 1.1733 and the distance to reson ance is A = 
PclPb Gl7 — 1 = 0.0057, using periods listed bv ICarter et al.l 
(l2012ll . If Kepler 36b was once located 0.4% further away 
from the star, we find A = 0.0001, which is much deeper 
inside the resonance. However, the natural outcome of the 
torque we have derived here is to move the planet outward, 
not inward. In principle, inward migration could occur In 
the case of a retrograde zonal wind. Either, a wind from 
Kepler 36b did not exert a significant torque on its orbit, or 


the wind moved the planet only slightly inward rather than 
outward. 

Dissipation in planets from variations in the tidal 
force from the central star can slowly cause orbits 


of planet pairs in mean motion resonance to diverge 
IPa^aloizm^ 20n Ddisle_et^ 20121: Lithwick &: Wull2012 : 


iBatygin fc Morbidellil 20131 : Lee et al. 2013h . Tidal dissi¬ 
pation is a mechanism that could account for the many 
near resonant Kepler planet pairs (ILithwick fc Wu 20li 
iBal ^in fc Morbidellil [201^ Leeet^^ay" 2013h . Quillen et alJ 
I 2013li found (following iLee et al.ll2013l ~) that it was difficult 
to account for the distance of the Kepler 36 system from the 
7:6 resonance (assuming it was originally left there) using 
tidal forces alone (though the effect of larger planet radii 
and masses at earlier times and prior to evaporation was 
not taken into account). 

Some models of tidal evolution p redict low eccentricity 
after escape from resonance (see, e.g.. ILithwick fc Wu 20^, 
with the exception of the scenario by Delisle et al.l 201^ . 
However the Kepler 36 system is near some second order 
resonances and, if there was divergent evolution a resonance 
crossing could perhaps have lifted the eccentric ity of the two 
plane ts to near their current values ~ 0.02 llCarter et al.l 

l2012h . 


While second order resonance crossing (due to drift 
caused by an evaporative wind) presents a possible explana¬ 
tion for Kepler 36’s eccentricity, lower index resonances such 
as the 3:2 are more distant from second order resonances. 
Only if the wind increases the eccentricity could non-zero 
eccentricities in lower index resonances be related to orbital 
variations associated with an evaporative wind. 


6 SUMMARY AND DISCUSSION 

In this paper we have studied the torque exerted by an evap¬ 
orative wind on planets with properties of super-Earth or 
hot Neptunes. 

Evaporative winds driven from close-in planet by EUV 
and X-ray radiation from a star are unlikely to be isotropic. 
If the upper atmosphere of a wind undergoing heavy mass 
loss is super-rotating, then we estimate that there would be a 
lag on the dawn side between the time the upper atmosphere 
is illuminated by starlight and when a newly launched wind 
becomes transonic. This lag results in a torque on the planet. 
A planet that loses a significant fraction of its mass could 
drift by a small amount in semi-major axis. By estimating 
wind launch timescales we find that only in a narrow regime 
of planet radius, mass and UV -I- X-ray flux can the lag 
angle be significant. Consequently only in rare cases could 
a planet’s orbit (or eccentricity) be affected by the wind’s 
torque. 

Our ID, time-dependent, hydrodynamical simulations 
have shown that for a close-in, Neptune-like planet, the 
anisotropic torque exerted by the wind does not strongly 
affect the semi-major axis of the planet. For the cases we 
have considered here, a Neptune-like planet on a short pe¬ 
riod orbit could undergo a semi-major axis drift of the order 
of 1%. However, this small effect has slightly more significant 
consequences in the case of Kepler 36b. The small change 
in semi-major axis that the wind could cause might be im¬ 
portant in sculpting the near-resonant configuration of the 
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Figure 6. 2D map of sonic radius divided by planetary radius for a pre-evaporation Kepler 36b planet, as a function of longitudes 
(day-side and night-side, 0 is dawn, tt is dusk) and latitudes (from south pole to north pole, 0 is at equator). 



r [planetary radius] 

Figure 7. Lag angle 6 as a function of the radial distance r from 
the planetary surface, for a pre-evaporation Kepler 36b planet. 
The wind is rapidly launched in the day-side, and its density 
decreases at large radii. On the night side, there exists a quasi 
steady-state and isothermal wind, with a significant density at 
larger radii. This could explain the small, and sometimes negative, 
lag angle. 


system. We find that the amount of drifting in semi-major 
axis caused by an anisotropic wind toque could possibly ac¬ 
count for Kepler 36 b,c’s distance from the 7:6 resonance. 

More generally, this mechanism has important conse¬ 
quences for close-in resonant pairs of planets. If the wind 


circulates in the prograde direction, the torque arising from 
the evaporative wind would cause outward migration of the 
innermost planet. Assuming that the outer planet is less af¬ 
fected (as it receives a weaker flux) then the pair would be 
pushed deeper inside resonance. Therefore, any mechanism 
that predicts resonance repulsion should be strong enough 
to counteract the torque arising from an evaporative wind. 
If the wind circulation was retrograde, then the innermost 
planet would be pushed inward. We have found that for the 
case of Kepler 36b, it is enough to bring it to its current 
location just outside the 7:6 mean motion resonance. Hence, 
there is strong motivation to do global circulation models 
for strongly irradiated hot Neptunes. We note that for low 
order resonances, such as the 2:1 resonance, the effect of the 
torque on A (distance to resonance) might be weak because 
of the large orbital separation between the two planets (com¬ 
pared to 7:6 resonance, for instance). However, we point out 
that M type stars have a strong X-ray flux for about 1 Gyr, 
so the time-integrated effect for planets around such stars 
will be more important. 

The values we have estimated for the lag angle depend 
on the zonal speed we use (here 1 km/s). If the zonal winds 
were rotating slower, the lag angle would also be smaller, 
which would minimize the influence of the torque. In this 
sense, the lag angle that we derive here can regarded as an 
upper limit on how an evaporative wind can affect a planet’s 
orbit. Global 3D simulations of evaporative atmospheres 
could also put more accurate constraints, since our 1-l-lD 
simulations can only be considered illustrative of the kind 
of anisotropic winds that will arise. 3D simulations would 
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not suffer from the uncertainly in the lag angle that arises 
from our approach. 
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APPENDIX A: ORBITAL EVOLUTION 

In this section we compute the rate of change of semi-major 
axis and eccen tricity c aused by a wind’s momentum loss 
rate. Following iBurns IGizi), we average the force per unit 
mass exerted by the wind on the planet over the planet’s 
orbit. 


A1 Semi-major axis evolution 

A wind that carries momentum at a rate P exerts a force 
per un it mass Fu, = P/Mp on the planet. Following iBurnsI 
lll976ll we decompose the force per unit mass from the wind 
into components 

F^=Te + Rr + Nz (Al) 

in spherical coordinates defined by the planet’s position rf 
and with z aligned with the orbital angular momentum vec¬ 
tor. Here T, R, N have units of acceleration. The work per 
unit mass (change in orbital energy per unit mass) is 

E„ = vF^=rR + rOT. (A2) 


The change in orbital energy directly gives an estimate for 
the drift rate in semi-major axis 

2a^ 

= + (A3) 

At low eccentricity, the timescale for drift in semi-major axis 
is 


Ta = 


iGM* 


2ra2 


(A4) 


where n is the planet’s mean motion. Because hp is non-zero 
at zero eccentricity, we ignore the sensitivity of the migration 
rate to eccentricity. 


A2 Eccentricity evolution 


For the eccentricity evolution 


^ dt 2 


2HEo H 


H + 


n 

2HEorT 


-Eo 


tt2 

+ - (r-R + reT) 

n 


where the angular momentum per unit mass reads 


(A5) 


H = 


= napyz. 


(A6) 


Here e is the planet’s eccentricity and y = Vl — e^. Whereas 
hp contains a non-zero eccentricity independent term, eccen¬ 
tricity damping or growth only takes place if there is eccen¬ 
tricity and both terms in eauation lA5l are needed to calculate 
this rate. The torque H = r x F, and the magnitude H = rT 
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in terms of the tangential force component. The radius from 
the star is 


r = ap(l — ecosE) 


(A7) 


with E the eccentric anomaly. The mean anomaly M and 
the eccentric anomaly are related by Kepler’s equation 

M = E — esinE. (A8) 

The mean motion 

n = E(1 — ecos E). (A9) 


Inserting these two equations into our expression for radius, 
we get 


. apesmEn 

r = apCSinEE = 

1 — e cos E 


Using H = r^6 we find 


rO = 


naly 


napy 


1 — e cos E 


(AlO) 


(All) 


Putting expressions that depend on the eccentric 
anomaly into equation (IA5II . we find for the eccentricity evo¬ 
lution 


-1 

n 


dt 2 



GM, 


—r(l — e cos E)y + 


Ty^ 

1 — e cos E 


+ 


Re sin Ey^ 
1 — e cos E 


(A12) 


This is equivalent to equation (28) bv iBurns IdTaiS) but this 
is written in terms of the eccentric anomaly so it can more 
easily be integrated. 

Let us assume that R and T both depend on radius in 
the same way so that 


T = R sin 


(A13) 


with 0U, the angle of the integrated outflow momentum. This 
means that T < R. 


_i d 


Raz 


dt 2 GM, 


sinO^y ( —(1 — ecosE) + 


V 


1 — e cos E 


+ 


e sin Ey^ 


(A14) 


1 — e cos E 

Differentiating Kepler’s equation 

dM — ndt = dE{l — ecosE) (A15) 

and we average equation (IA14II over a rotation period 


[sin 6njy{2e cos E — e^ cos^ E — e^) 

4-p F?/^l 


There is no eccentricity damping unless there is some eccen¬ 
tricity. Hereafter we drop the brackets and discuss the orbit 
averaged eccentricity evolution. 

Assume that with a low order eccentricity expansion we 
can write 


R{E)R{1 + Tc cos E + rs sin E) (A17) 
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where and Ts are unit-less but could depend on eccentric¬ 
ity. If the momentum flux R{E) oc r~^ then to first order 
R{E) oc (1 — Pecos E) and rs = 0. If there is a seasonal lag, 
then we might get non-zero Vs term. Integrating equation 
(IA16II we find 


-1 

n 


d 

dt 2 


Rap 

gm/ 




(A18) 


and we have dropped terms higher than second order in 
eccentricity because equation (IA17I) is a low order approxi¬ 
mation. 

If the wind momentum flux depends inversely on radius 
then Tc ~ e and 


n 


-1 


d 

dt 2 


RUp 

GM^ T 


sin 6w 


(A19) 


and we expect eccentricity damping. If the wind momentum 
flux depends inversely on the square of radius then ~ 2e 
and 


_i d 

dt 2 


Ral 

GM, T 


sin 9ui 


(A20) 


and w e expect an increase in eccentricity. I Owen fc JacksonI 
l|2012tl predict that X-ray driven winds have mass loss rate 
proportional to and EUV driven winds have mass loss 
rates proportional to r~^. This suggests that for planets 
within 0.1 AU of a star, the wind could cause the orbital 
eccentricity to increase whereas external to this the eccen¬ 
tricity would be damped. 

And in both cases the timescale for eccentricity varia¬ 
tion is similar to 


de'^/dt 


1 GM, 
nsinS™ Ra^ 


(A21) 


Our expression for the eccentricity damping (or in¬ 
crease) timescale and eccentricity evolution rate depend on 
Ra^/[GM,) = R and R is in units of acceleration as it is is 
force per unit mass. This implies that 


5 Mvw 

Rop v-w M 

GMt Vc MpU 


(A22) 


(A23) 


which depends on the ratio of the wind outflow velocity and 
orbital speed and the mass loss rate related to the fraction 
of mass lost per orbital period. The torque we previously 
estimated is = MpT, giving 


sin0^ = . ™ (A24) 

Mu™ 


Putting these together , we find: 

_ _ 1 Uc Mp 

de^/dt sin0™ u™ M 


(A25) 


and Tf. ^ Ta- We find that the rati o commonly used to de¬ 
scribe migration rates, K = Ta/re l|Lee fc Pealell20fl3 l. is of 
order unity for evaporative winds. 
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